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Abstract 

This paper is concerned with a lattice model which is suited to square-rectangle 
transformations characterized by two strain components. The microscopic model in¬ 
volves nonlinear and competing interactions, which play a key role in the stability 
of soliton solutions and emerge from interactions as a function of particle pairs and 
noncentral type or bending forces. Special attention is devoted to the continuum 
approximation of the two-dimensional discrete system with the view of including 
the leading discreteness effects at the continuum description. The long time evolu¬ 
tion of the localized structures is governed by an asymptotic integrable equation of 
the Kadomtsev-Petviashvili I type which allows the explicit construction of moving 
multi-solitons on the lattice. Numerical simulation performed at the discrete system 
investigate the stability and dynamics of multi-soliton in the lattice space. 


1 Introduction 

A lot of interest has recently been devoted to spatio-temporal patterns as well as the asso¬ 
ciated defects and dynamics such as standing-wave patterns, localized structures including 
solitons or oscillating patterns. These structures become fundamental in the study of phase 
transitions which are usually accompanied by the appearance of defects: dislocation mo¬ 
tions, grain boundaries, domain wall structures and twinnings [l]-[3]. One of the aims 
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of the present research is to understand how solitons arising at the microscale (ie at the 
level of the lattice model), are able to organize the system at the macroscale and what 
the selection properties of the nonlinear structures are. The present work is particularly 
motivated by the existence of solitonic structures occurring in phase transformations in 
crystalline alloys. 

The same lattice model has been successfully used to examine the formation of localized 
strain patterns decaying in all directions [4]-[6]. Then the corresponding research results 
in the partial softening of the transverse-acoustic phonon branch at a nonzero wave num¬ 
ber; the positive curvature of the dispersion branch at the long wavelength limit; and the 
shearing motion of the atomic planes along the stacking direction leading to spatially ar¬ 
ranged structures made of martensitic twin bands and the existence of strain solitary waves 
describing the coherent movement of martensitic domains. More precisely, the nonlinear 
dynamics of the two-dimensional model allows us to examine the stability of the lattice 
using numerical simulations which have shown the formation of localized strain structures 
emerging from an instability mechanism [4, 6]. 

In the present paper we continue the investigation of the properties of the two-dimensional 
lattice model in order to study the existence, stability and dynamics of lattice multi- 
soliton conhgurations. In fact, the study provides the critical values of the coefficients 
of the lattice model for which lattice solitons exist and move along the plane, using the 
fact that the continuum limit of the long time evolution of the patterns is governed by an 
asymptotic equation of the Kadomstev-Petsviashili I type. Recently in [7], the Kadomstev- 
Petsviashili equation has been derived for quasiplane waves by considering oscillations of a 
two-dimensional square lattice array of atoms with interactions only between the nearest 
neighbors. 

This paper is organized as follows: in section 2 we introduce the lattice model, while 
in section 3 we deal with its continuum approximation and show that it leads to the 
Kadomstev-Petsviashili I equation for long time evolution. Finally, in section 4 numerical 
simulations of the discrete model based on the Kadomstev-Petsviashili I multi-solitons are 
performed and their dynamics are discussed in some detail. 

2 Lattice Model 

Let us consider an atomic plane made of squares parallel to the i and j directions presented 
in Figure 1. Such a lattice model can be extracted from the cubic lattice of crystalline 
alloys (like the fee symmetry of In-Tl, Fe-Pd and other crystals). The model describes a 
cubic tetragonal transformation of the lattice. A particle of the lattice plane is located by 
(f,j). After deformation of the lattice, the particle undergoes a displacement dehned by 
Uij = u{i,j) along the i and j directions. 

The particles interact via two kinds of interatomic potentials: (i) interactions between 
hrst-nearest neighbors considered as functions of the particle pairs in the i and j diago¬ 
nal directions and (ii) interactions involving noncentral forces or three body interactions 
between the hrst-nearest and second-nearest neighbors in the i and j directions. The 
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Figure 1: Two-dimensional lattice model with the detail of interatomic interactions between 
the hrst and second neighbors in the i,j direction. 


potential describing the first interactions possesses stable, unstable or metastable states 
accordingly to a control parameter which is connected with temperature [8]. The latter 
interactions amount to describing, at the microscopic level, the resistance of the crystalline 
cell to twisting and bending [4, 5] and thus provide competing interactions [9, 10]. The 
noncentral interactions are of particular interest for the competing interactions and sta¬ 
bility of the nonlinear structures. On using the invariance of the lattice energy under 
translations and rotations, the following potential functional was introduced in [6] 




(id) ' 




( 2 . 1 ) 


where the discrete deformations are dehned as 


^i,j ^id 


^i,j ^id 


( 2 . 2 ) 


and the potential $ is given by 

^ ^ a2 Sl^. (2.3) 

For simplicity, the lattice energy (2.1) has been chosen to be dimensionless. The first 
and second terms in (2.1) represent the nonlinear and linear potentials coming from the 
particle pair interactions where (ai, 0 : 2 ) and (3 are the lattice force coefficients for the 
longitudinal and shear deformations, respectively. The third and fourth parts of (2.1) hold 
for the actions of the noncentral interactions in the i and j directions. The interactions 
are characterized by the parameters <5 and rj for the actions between first- and second- 
nearest particles, respectively. On the other hand, the operators A]^ and A^ hold for 
the forward hrst-order hnite difference in the i and j directions: = 5'j+ij — S'jj 

and A^Gjj = Gjj+i — Gjj. Remark: Due to (2.2) the noncentral interactions in (2.1) 
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are of the form Mj-ij + Mi+ij — for the first-nearest neighbor interaction and of the 
form Ui- 2 ,j + for the second-nearest neighbor interaction (same for the Gij 

deformation). 

In addition, a one-dimensional version of the model can be obtained, reduced from the 
complete two-dimensional one, in which the shearing motion of the atomic planes along 
the stacking direction is modeled by arrays of martensitic and austenitic solitary waves. 
This reduced system has been examined in detail, and it has been observed that localized 
structures emerge of arrays of elastic solitary waves [4, 5]. 

By introducing the kinetic energy associated with the displacement Uij (for unit mass) 
as 

(2-4) 

hj) 

the corresponding difference-differential equations of motion for Uij, deduced from the 
Hamiltonian H = T + V, are given by 

Sr,, (2.5) 

where the discrete stresses are dehned as 

= cTi^j - Al ( 2 . 6 ) 

= P Gi,j - XTij (2.7) 

(Tjj Oil ^i,j (^•^) 

XLij = ^i,j + V [^i+2,j + -|- QSij + ASi-ij + Si-2,j]} (2-9) 

XTij = ^i,j V [^i,j+2 ~‘t~ 4:Gij+i + 6Gij + AGij-i + Gij-2\} ■ (2.10) 

Equations (2.6-2.7) correspond to the discrete macroscopic stresses due to the fact that the 
deformations Sij and Gij are functions of the discrete displacements Uij. Note that the 
stress (2.8) which follows from the potential <h since aij = d^/dSij is nonlinear in terms 
of the strain Sij. Also, the microscopic stresses (2.9-2.10) emerging from the noncentral 
forces are functions of the discrete variations of the deformations Sij and Gij in the i and 
j direction, respectively. 

Due to the strongly nonlinear nature of the problem, these equations are not managable 
except for the linear problem which has been examined in Ref. [11]. In this paper, these 
equations are solved using numerical simulations with appropriate initial and boundary 
conditions. More precisely, using the continuum approximation a quasi-continuum model 
has been obtained which includes the leading discreteness effects and allows us to investi¬ 
gate the existence of soliton configurations and study their dynamics. 

3 Continuum Approximation 

In order to describe the lattice dynamics at the quasi-continuum scale we assume that the 
discrete functions are slowly varying over the lattice spacing. In fact, both the deformations 
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(2.2) have been expanded using Taylor’s series up to third order, that is S' = hux — ^Uxx + 
^Uxxx (where h is the lattice spacing); while all other terms have been expanded for long 
wavelength up to fourth order. After some classical algebras the continuum equations of 
motion for the deformation u{x,y,t) are 


^tt ^xx ^2 (^ 'jxx T ^xxxx T ^yyyy A f^^yy (3.11) 

for dimensionless space unit, ie rescaling x ^ x/h and y ^ y/h. The coefficients 5l and 
5t are given in terms of the coefficients of the model, ie 

= ^ - 5 - 16r/, 5T = ^-5-lQy. (3.12) 

3.1 Asymptotic Model 

In order to understand the evolution of the localized structures over a large scale of time 
of the order of and for a weakly nonlinear medium, we consider the asymptotic equa¬ 
tion derived from the continuum equation (3.11). By assuming that the contribution of 
the nonlinear term is weak throughout we rescale the nonlinear and dispersive terms by 
introducing a small parameter e < 1 as follows 

Oi2 = £C(, — £^t- (3.13) 

Then for large times (of order £~^), the asymptotic expansion of the displacement held is 

u{,x, y, t) = uo{^, Y, t) + £Ui{x, y, t) + O(e^) (3.14) 

where ^ = x — ct is the stretch phase variable, Y = £^^‘^y is the stretch transverse variable 
and T = £t is the slow time variable. 

Using equations (3.13) and (3.14) and keeping terms of order e only, equation (3.11) 
becomes 

2cMo,e “ " (“U? + “oyy = Ml,, - (3.15) 

where we have set the sound velocity in the x direction to be equal to c = ^Jol{ (for ai > 0). 

Then the secularity condition [12] implies that both the left- and right-hand side of 
(3.15) have to equal zero. Thus, the right-hand side leads to the standard two-dimensional 
linear wave equation for ui, while the left-hand side gives the long time evolution of Uq 
dehned by 

(moj, 4- a{ul)^ + = Puoyy- (3.16) 

Note that the following change of variables and of parameters have been considered T = 
—r/2c and Sl = —6l in order that (3.16) transforms into a standard form, and that for 
specihc choices of the parameters d, Sl and P > 0 equation (3.16) becomes the Kadomtsev- 
Petviashvili I equation [13] (see later). Note that equation (3.16) was hrst derived in [6] 
using Fourier images. 
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3.2 The Kadomstev-Petsviashili Equations 

Nonlinear quasi-one dimensional waves (with y much larger that x) in a weakly dispersive 
medium are described by the dimensionless Kadomstev-Petsviashili (KP) equations 

3('U tixxx^ ^ ^yy 


where the ± on the right-hand side of (3.17) is determined by the dispersive property of the 
system; ie the upper sign is usually refered to as positive dispersion and the corresponding 
equation is known as KPI (which is the case we study here). 

The exact multi-soliton solutions of KPI can be constructed by different methods (see, 
for example, Refs. [14, 15]). For our purposes it is convenient to write them in the Hirota 
form, ie 


u{x,y,t) = 2 


dx‘^ 


In I 


(3.18) 


where 0 is the determinant of a 2 n x 2n matrix (for n number of solitons) given by 


0 = det 


{x -h PkV + pit + Ok) 6ki + (1 - hi) 


2V3i 

Pk-Pi 


(3.19) 


Here 5ki is the Kronecker symbol, the indices take values k,l = 1,2, •••,2n, while the 
complex constants pk and Ok determine the velocity and the phase of each soliton, where 
Pk+n Pk and Ok-\-n Ok- 

The one soliton solution given by (3.18-3.19), when pi = i and 0i = 60, has the simple 
expression 


u{x,y,t) = 4 


-{x-t + 60)^ + y^ + S 
{x — t + 60)^ + y"^ + 3 


(3.20) 


The corresponding conhguration consists of a soliton which for t = 0 is situated at the 
points {x, y) = (—60, 0 ) while for t 7 ^ 0 travels along the x-axes without changing its shape 
and with constant velocities: {vx,Vy) = (|piP, —2 3fJ(pi)). 

In addition, the corresponding multi-soliton solutions (3.18-3.19) are asymptotically 
free at t ^ ±cx) and move along their unperturbed trajectories without changing their 
shapes and initial parameters after collision, ie they scatter trivially. Figure 2 represents 
snapshots of the two soliton solution (3.18-3.19) for pi = i, P 2 = 2i and 0i = O 2 = 0. 

Note that in the framework of plasma physics the classical KP equations suffer trans¬ 
verse instabilities accordingly to the sign of the dispersive terms as well as possessing soliton 
solutions [16]. 


4 Numerical Simulations 

In this section, we investigate numerically the existence of lattice multi-solitons using the 
fact that the lattice model asymptotically leads to an integrable equation. We use a 
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Figure 2: The two soliton solution (3.18-3.19) of the KPI equation at different times with 
trivial scattering behaviour. The taller soliton passes through the smaller one with no 
phase shift or radiation. 


numerical scheme by directly considering the lattice equations for the displacement Ui^ 
given by 

(5 ('U^-|-2,j “1“ QUi^j Ij' ~\~ Ui—2J^ ^ (^^,7+2 “t” GUi^j AUi j—\ Ui^j—2^ 

Tj ('U2-|-4j' 4'U2-|-2,j T d'U^_2,j “t“ Tj ('U^j-|_4 4'U^j-|_2 “I” 2 T 

(4.21) 

which follow from (2.2) and (2.5). The numerical simulations are performed by employing 
a Runge-Kutta method of fourth order and by imposing periodic boundary conditions, ie 

'^i+N,j '^i,j+N ^i,j (4.22) 

where N is the number of particles along the boundaries of the square. 

Note that it is not obvious if the solutions of the KPI equation will also be a solution 
of the lattice model. However, it was shown in [17] that this approach has been success¬ 
fully considered in a one-dimensional lattice model. Specihcally, when a one-dimensional 
lattice model with long-range interactions was considered (which, in the continuum, keeps 
its nonlocal behaviour) the long time evolution of the localized waves is governed by an 
asymptotic integrable equation the so-called Benjamin-Ono equation which allows the ex¬ 
plicit construction of moving kinks on the lattice. Accordingly, in this section we investigate 
the dynamical behaviour of the lattice model (4.21) numerically, using as initial conditions 
the one and two soliton solution of the KPI integrable equation given in section 3.2. 

By comparing the coefficients and the variables of equations (3.16) and (3.17), it is easy 
to see that, (3.16) transforms to (3.17) when 

(4 = 3, 

(3 = 1, 

Sl = 1 , 


X ^ ^ = X — ct 

y^Y = 

t^T = -etl2c. (4.23) 
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Figure 3: The one soliton solution (4.25) for e = 0.1. 


Then it is a matter of simple algebra to observe that, due to (4.23), the parameters of the 
discrete model given by (3.12-3.13) are equal to 

02 = 3e, 5 = — -16r] - e (4.24) 

and, thus, the only arbitrary parameters of the lattice model are the following three: e, 
c and Tj. Due to numerical simplicity, in all our simulations we have chosen the following 
hxed values for: (i) the small parameter related with the weak nonlinearity e = 0 . 1 , (ii) 
the lattice spacing h = 1, (iii) the total number of the lattice points N = 100, (vi) the 
time step dt = 0.1 and (vi) the total number of time steps equal to 5000 points. However, 
different values for the above parameters give qualitatively the same results as long as e < 1 . 

• Dynamics of One Soliton 


Firstly, we investigate the time evolution of a single solution where the conhguration 
given by (3.20) have been considered under the change of variables (4.23). More precisely, 
the initial conditions for the displacements and the velocities of each lattice particle, in the 
case of the one soliton (3.20), are given by the analytic expression 


^ — {x — ct + et/2c -|- 60)^ -|- ey‘^ + 3 
[(x — ct + etj^c + 60)2 _j_ ^y2 _|_ 3j2 


(4.25) 


and its time derivative, respectively. Figure 3 illustrates snapshot of the displacement 
Uijit) at t = 0 . 

We have run our simulations for different values of the parameters (c, 77 ) and investigate 
whether the initial lattice soliton relaxes and propagates in the lattice space with small 
(or no) oscillations. This phenomenon occurs, only, for specihc range of values of the 
aforementioned parameters which give qualitatively the same results. In particular, due 
to [ 11 ], we know that the noncentral interaction parameter have to be negative and small. 
Our numerical study have shown that for values of the parameter 77 < 1/60 and when the 
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soliton velocity takes values between 1 < c < 4 the loss of energy by radiation of small 
amplitude waves is small. 

In figure 4, we display a full three-dimensional plot corresponding to the displacement 
Uij at four different times f = 0,14.5, 29.5, 50.5. The initial conditions were created from 
one soliton placed at the position xq = —60. Figure 4 shows that the lattice soliton is 
stable throughout the numerical simulations; more precisely, its size is constant as it moves 
towards the plane without emitting any radiation. This process continues through several 
cycles, due to the periodic boundary conditions, without the occurrence of instabilities 
since the conhguration settles to a single soliton which move along the lattice grid. 

Qualitatively, similar results have been found and studied in great detail in [17] for 
one-dimensional lattice kinks which are solutions of a nonlocal discrete model with long- 
range interactions. In that case, the lattice radiation was inversly proportional to the kink 
thickness. 

• Two Soliton Scattering 

Next, we discuss the results of a numerical evolution of the full time-dependent lattice 
equations (4.21), in order to investigate the interaction and scattering of two solitons. 

We take two solitons initially located at the origin (represented in Figure 2) and evolve 
the equations of motion for different values of the soliton velocity c and of the noncentral 
force parameter rj. For values (approximately) between the following range: 4 < c < 8 
and —1/30 < rj < —1/360 the numerical simulations give qualitatively the same results, 
represented in Figure 5. 

From Figure 5, we see that the two solitons which initially (at t = 0) form a single 
large soliton travel with opposite velocities away from each other and then, due to the 
periodic boundary conditions, reappear at the edge of the lattice. Now, the two solitons 
head towards each other until they merge to form a single soliton again (at t = 48.5) while 
after the collision, they continue their path with no phase shift or emission of radiation, 
ie they scatter trivially as the KPI solitons. In fact a small amount of energy is radiated 
during the soliton motion as a result of the highly nonlinear deformation the solitons suffer 
in adjusting to the lattice grid. This radiation can be dealted with numerically by applying 
absorbing boundary conditions at the edge of the grid. The radiation emission is eliminated 
as the soliton velocity c increases. This scattering behaviour is continuous throughout for 
a long time; specihcally, we have run our simulations up to t/ = 500 and although the 
lattice grid is covered from radiation the two solitons are easily distinguished. 

To conclude, this approach can be extended to higher soliton numbers and we expect 
the results to be qualitatively like the ones obtained in the hrst two cases of the one and 
two soliton dynamics. 
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5 Conlusions 


Tha main objective of the present work is to examine the long-time evolution of nonlin¬ 
ear localised waves of soliton type propagating on a two-dimensional lattice model. The 
starting lattice model involves nonlinear (nonconvex lattice potential) and competing in¬ 
teractions which follow from noncentral interactions and play a key role in the existence 
and stability of the localised waves. Physically, the model describes the mechanism of 
micro-twinning of small ferroelastic domains in alloys suffering phase transformations like 
ferroelastic materials, martensitic transformations in shape memory alloys, etc. 

The most interesting point of the study is that, on the basis of a multi-scale technique, 
an asymptotic model for the long-time evolution of the nonlinear waves has been derived. 
This model is then governed by Kadomtsev-Pietvishviali I equation which possesses soliton 
solutions. In order to check the analytical conjecture given by the asymptotic model, some 
numerical simulations have been performed (directly) on the discrete system. In fact, the 
numerical results describe the time-evolution of one- and two-soliton solutions. It appears 
that the soliton solutions thus generated are particular robust and stable which implies that 
the prediction provided by the asymptotic model is a good approximation of the solution 
of the discrete system. Therefore, we conclude with the observation that the KPI equation 
is an efficient asymptotic model to predict the existence of localised objects or nonlinear 
excitations of soliton type in a rather complex physical system. 

Further extensions of the present work can be investigated. In particular, the dynamics 
of localised waves for a lattice model including higher order nonlinear interatomic terms in 
the potential (fourth order nonlinear term in discrete deformation, see equation (2.3)) can 
be considered; while it is worth investigating the influence of applied forces and damping 
on the nonlinear wave propagation and stability. 
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